Vortex String Formation in a 3D U(l) Temperature Quench 
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We report the first large scale numerical study of the dynamics of the second order phase transition 
of a (7(1) \(j} 4 theory in three spatial dimensions. The transition is induced by a time-dependent 
temperature drop in the heat bath to which the fields are coupled. We present a detailed account of 
the dynamics of the fields and vortex string formation as a function of the quench rate. The results 
are found in good agreement to the theory of defect formation proposed by Kibble and Zurek. 
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Topological defects are of fundamental importance to 
the (thermo) dynamics of phase transitions in some of the 
most fascinating materials in the laboratory - eg. super- 
fluids, type II superconductors, liquid crystals - and pre- 
sumably also to symmetry breaking phase transitions in 
the early universe Q . When present at low energies they 
constitute the last traces of disorder inherited from high 
temperatures and thus determine the non-equilibrium 
evolution of many systems. 

The theory of defect formation combines the realiza- 
tion, due to Kibble of the paramount role of causality, 
with the calculation, due to one of us ||, of the actual 
size of the causally independent domains in a second or- 
der phase transition. Experiments in liquid crystals 
and in 3 He || lend support to the resulting theory of 
the dynamics of second order phase transitions. The evi- 
dence from 4 iJe experiments |^| is more ambiguous while 
superconductor and Bose-Einstein condensates may offer 
exciting future possibilities 0]. 

Laboratory experiments were, so far unable to test the 
key theoretical prediction - the scaling of the initial den- 
sity of defects with the rate at which the phase transi- 
tion takes place. Numerical studies carried out to date 
were limited to 1 and 2D systems, where they have con- 
firmed the scalings predicted by the theory p. More- 
over, they have focused on transitions induced by the 
explicit change of the mass term in the potential which 
governs dynamics of the order parameter field (0), while 
the temperature T of the heat bath to which <f> is cou- 
pled was kept constant and relatively small. This sim- 
plification is not necessarily unrealistic - phase transi- 
tions in 4 He can be induced without significant changes 
in temperature, purely by changing pressure. However 
such 'pressure quenches' are far less common than tem- 
perature quenches, and temperatures are rarely small. 
Moreover, the analytic approach to the pressure quench 
problem seems to be within reach, although only in the 
limit of very small temperatures. Even in this limiting 
case unanimous conclusion has not been reached [p|-pd]|, 
further obviating the need to study quenches numerically. 



Indeed the emerging understanding of the defect pro- 
duction could be called into question on the grounds that 
the present estimates completely ignore the inevitable 
Ginzburg regime. There, just below T c , fluctuations 
can rearrange large spatial regions, which could destroy 
proto-defects produced by freeze-out 0, or create them 
at densities set by the Ginzburg length, as it was origi- 
nally suggested 0] and is still occasionally argued || . 

In this letter we perform the first large-scale numerical 
studies of a temperature quench in 3D U(l) symmetric 
A<fi . This is the familiar Ginzburg-Landau model de- 
scribing the free energy of a neutral system, (eg. 4 He), 
which as is well known displays a true second order phase 
transition in 3D, with the establishment of long range 
order at low temperatures. This implies in particular, 
contrary to the case of lower dimensions, that no defects 
(vortex strings) can exist in equilibrium at sufficiently 
low temperatures in any causally connected volume [fl2| . 

In order to implement the temperature quench we 
evolve the fields according to: 



(d 2 - V 2 ) & - m 2 ^ + + <$) + 



(1) 



where i,j € {1,2} and i ^ j in Eq. (|l]). Fi(x) is the 
Gaussian noise characterized by 

(Ti(x)) = 0, (F^r.Or')) = 2 V T(t)S ij S(x - x'), (2) 

where x, x' denote space-time coordinates. We allow the 
fields to thermalize above the transition and proceed to 
quench the system by changing the noise temperature as 



T{t) =T C -T a 



t 



(3) 



The time tQ controls the rate of the quench. For t > 

In the numerical evolution we take the grid spacing 
Ax = 0.5 and At = 0.02. All results shown are for com- 
putational domains of size N 3 , with N = 128 — 160 and 
rj = 1. Strings are detected by integer windings of the 
field phases around lattice plaquettes. Vortices hence 
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found are connected by enforcing flux conservation on 
each unit volume ]T^ | . As a result of boundary conditions 
all strings are closed. Above T c , and immediately below 
strings are little more than non-perturbative field fluc- 
tuations. Below T c they gradually acquire stability and 
can be regarded as either vortex lines or cosmic strings. 

Another useful quantity is the kinetic temperature, 
Tk, defined as the average kinetic energy per degree of 
freedom. In equilibrium the canonical momentum distri- 
bution is purely Gaussian and {\ir{x)\ 2 )dP x — IT. We 
generalize this for situations away from equilibrium as 



(\7r(x 1 t)\ 2 )d D x = 2T K (t). 



(4) 



All temperature quenches are started at Tq = 1.91T C . 
Even though the choice of initial (high) temperature is 
somewhat arbitrary, it is important that it is sufficiently 
high that the length density in long strings in each com- 
putational domain is substantial |l2|. 

In the immediate vicinity of T c for a system undergoing 
a second order phase transition, the dynamics of <f> are 
subject to critical slowing down. This leads to the esti- 
mate of expected density of defects through an argument 
||, which we briefly reproduce below. For the dynamics 
of Eqs. in the over-damped regime where the first 

time derivative dominates, the characteristic time scale 
t over which the order parameter can react is given by 



Tl 



V 



(5) 



where v and z are universal critical exponents, and the 
relative temperature e — ^ — 1 = = ^y 1 - The 
quench timcscale tq is a 'natural units' rescaling of tQ, 
T Q = 5§*Q- 

This over-damped scenario, valid when t] 3 tq > 1 Q, 
is presumably more relevant for condensed matter ap- 
plications and will the the focus of the present Letter. 
Cosmological order parameters may in contrast be under- 
damped -or in reality red-shifted jlj] , corresponding to a 
different dynamics than that of Eq. ([j]j^). 

The characteristic time scale of variation of e is 4 = t. 
We expect the system to be able to re-adjust to the new 
equilibrium as long as the relaxation time is smaller than 
t. Hence, outside the time interval [— £, t ] defined by 
the equation r(e(£)) = t, the evolution of <f> is approxi- 
mately adiabatic, and physical quantities associated with 
large length scales will approximately follow their (criti- 
cal) equilibrium values. The time 
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marks the borders between adiabatic and impulse stages 
of evolution of <j). In particular the correlation length £o 
associated the the connected 2-point function above the 
transition will cease to increase as £ = £o 771 o/l e | 1 ' once 
the adiabatic-impulsc boundary at t' — —t is reached. 



We expect, then, that the characteristic length scale 
over which <f> is ordered already in the course of the transi- 
tion will be the correlation length at freeze-out £ = 1/e", 



& = ( — T Q 



(7) 



On length scales smaller than £ the effect of criticality 
on the field dynamics is negligible. The initial density of 
vortex lines is then expected to scale with tq as 
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(8) 



where / ~ 0(1) is a dimensionless factor, parameteriz- 
ing our ignorance about the exact relation between do- 



main size and defect densities, and a 



1v 



Lattice 



measurements and renormalization group analysis yield 
v = 0.6705 and z — 2.03 for the very over-damped case 
when the second time derivative in Eqs. (0,§) is com- 
pletely negligible. With r\ = 1 this is not truly the case. 
We expect 1 < z < 2.03, implying 0.568 < a < 0.8, 
where the lower and upper limit refers to the under- 
damped case, respectively. This is different from mean- 
field exponents, v = 1/2 and z = 2, implying «mf = 0.5. 
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FIG. 1. a) The evolution of the total string density per 
plaquette n with Tk- ° denote thermal equilibrium densi- 
ties, b) vs. t. Solid lines denote quench trajectories for 
t Q = 4,8,16,32,64,128,256. 

In a temperature quench the system is evolved from an 
initial state at genuinely high temperature, with a bare 
negative mass squared. As a result the physical mass 
squared, m^, which takes into account the self-energy 
generated by the high temperature field fluctuations is 
positive and potentially large. The temperature quench 
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proceeds by the decrease of the external bath tempera- 
ture at a rate given by tq, according to Eq. (||). Ini- 
tially, the system locally re-thermalizes to the new lower 
temperature. Close to the critical point the physical 
mass squared approximately vanishes leading to the crit- 
ical slowing down of the field response over large spatial 
scales. This freezes the dynamics of the order parame- 
ter: It can no longer react to the systematic changes of 
thermodynamic or dynamical parameters, although slow 
drift under the combined influence of noise and damp- 
ing continues unabated, even in the large scale structure, 
including long strings. Critical slowing down has little 
effect over the small scale dynamics (fc 2 >> m^), which 
accompanies the externally imposed change of bath tem- 
perature. In this manner critical slowing down sows 
the seeds for the out-of equilibrium dynamics to follow. 
Fig. 0a shows the equilibrium string densities as a func- 
tion of Tff, and quench trajectories for several values of 
tq . In 3D Tk is dominated by small scales and generally 
follows faithfully the external bath T variation. 



scale structure in long strings and of small loops. For 
faster quenches this regime persists to much lower Tk- 
Fig. |l|a clearly shows the hierarchy of string densities for 
different tq's, at each given Tk- 

Fig. [l]b, shows the evolution of the order parameter 
amplitude (\<t>\)- The persistent cooling of small scale 
fluctuations eventually leads to a negative r/ip h which in 
turn triggers instabilities in the long wavelength modes. 
As these modes grow (quasi-exponentially) the U(l) sym- 
metry is spontaneously broken, coarsening the original 
field configuration. These instabilities are reminiscent of 
those in 'pressure' quenches [p|-{iT||. However here their 
generation through a negative mass squared is explicitly 
created by the delay in the fields response, rather than 
externally. The onset of instabilities defines in turn +t, 
which we read off from the minima of Fig. |b. Fig. |a 
shows the dependence of t on tq. We can independently 
confirm the theoretical scaling laws, Fig. ||b by examin- 
ing the tq dependence of e, computed as the ratio of Tk 
and T c at the minimum of {\4>\). 
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FIG. 2. a) The freeze-out time t as a function of the 
quench timescale tq, b) The relative kinetic temperature at t, 
1|. Slopes of the two lines are 0.445 ± 0.05 and 



0.565 ± 0.016, respectively, which compares favorably with 
theoretically predicted of 0.5. Points corresponding tq < 8 
exhibit saturation, and were ignored in the fitting. 

As predicted by the theory all quenched string den- 
sities follow the equilibrium trajectory at high temper- 
atures but start deviating just above the critical point, 
due to critical slowing down. The magnitude of the ef- 
fect is small and difficult to measure. After falling out of 
equilibrium the string densities do not freeze, but rather 
decay slowly while the temperature drops over a period of 
time. This decrease is due mostly to the decay of small 



FIG. 3. The string densities measured at (\cf>\) = 0.9, 0.925, 
0.95,0.975. The fits are to Eq. (|), for tq > 32, with 
a = 0.4296 ± 0.043, 0.4378 ± 0.0289, 0.5692 ±0.0256, 0.5600 
±0.0797 and = 11.11, 13.99, 12.29, 15.34, respectively. The 
average a = 0.4982 ± 0.079 and = 13.18 ± 1.78. 

Together the results of Fig. || confirm the theoretical 
scaling and determine vz — 0.82. We note however that 
data corresponding to the fastest quenches appears to 
asymptote to |e| = 1 (and, thus, cease following the pre- 
dicted power law). This is easy to understand, [||]: e, 
extracted from the data is the absolute value of the rel- 
ative temperature at the instant when the dynamics is 
re-started. But, e = Tk/T c — 1 by definition cannot fall 
to less than —1 (which happens when Tk = 0). Hence, 
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when in very fast quenches the evolution of <p re-starts 
only at t > tq, T c ^> Tk — T(tq) = 0, and the values of 
|e| pile up asymptotically near 1~, causing saturation. 

At t = t many small scale fluctuations still persist in 
the system obscuring the results in terms of string densi- 
ties. It is the subsequent out-of-equilibrium evolution of 
the fields, leading to spontaneous symmetry breaking and 
ordering, that reveals the string densities formed at the 
quench. While the large spatial scales (k 2 ~ 0) are unsta- 
ble and the corresponding modes grow towards \(j)\ = 1, 
the small spatial scales (k 2 » m 2 h ) are dissipated away. 
This leads to the emergence of the field configuration cre- 
ated by the critical dynamics on large spatial scales, and 
allows the surviving vortex strings to form, i.e. to acquire 
their low energy character of topological defects. These 
densities, as a function of tq, are shown in Fig. [|. As a 
criterion to the completion of the transition we measured 
n at t such that {\cj>\) = 0.9, 0.925, 0.95, 0.975. Apart from 
the cases of saturation, the density of strings formed fol- 
lows the theoretical predictions quite satisfactorily, for 
all choices of (|^|). The values of fi = O(10), are sim- 
ilar to 1 and 2D estimates ||, and may be sufficient to 
explain the non-appearance of vortex lines in the recent 
4 He experiment ff|. 



first. Nevertheless, the presence of noise and absence of 
the restoring dynamics implies that until t ~ +t a siz- 
able population of small loops can persist. Subsequently 
the system is dissipated further, the symmetry is spon- 
taneously broken and the fields order starting from the 
small scales. The result is that by the time (\(f>\) has 
come near its equilibrium low-temperature value, only 
long strings, imprinted by the critical dynamics on largest 
length scales, survive, stripped of most of their small scale 
structure. This evolution is shown in Fig. |]. 

Our numerical analysis lends strong support to the 
general picture of dynamical evolution of {\(j>\} in a sec- 
ond order phase transition, proposed some time ago [0, 
and partially confirmed by experiments j|] an d by simu- 
lation of 'pressure' quenches in low dimensional systems 
H . The key new aspects of this investigation are i) its 3D 
character, which has allowed us to ii) study a to data nu- 
merically unexplored temperature quench. In the present 
range of parameters the Ginzburg regime appeared to 
play no discernible role. We are currently investigating 
its effects on the decay of individual strings. 

We thank A. Gill, T. Kibble, P. Laguna, R. Rivers, and 
A. Yates for useful discussions. Numerical work was done 
on the T-division/CNLS Avalon Beowulf cluster, LANL 
and supported by DOE, contract W-7405-ENG-36. 
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FIG. 4. String length I distributions (n dl vs. 1) taken 
between +t and the 'time of formation' ({\4>\) = 0.95), for 
tq — 64. Data sets denoted by (*,A,D, •) correspond to 
increasingly later times. 

Finally, it is interesting to investigate what kind of 
strings are formed. At T above T c the 'ephemeral string' 
length distribution approaches a Brownian form [ jl2| . 
Long random-walk like strings of zeroes of 4> coexist in 
equilibrium with a sea of smaller strongly self-correlated 
loops. As the quench proceeds small scales are dissipated 
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